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ORIGINAL ARTICLE 

Schizophrenia-risk variant rs6994992 in the neuregulin-1 gene 
on brain developmental trajectories in typically developing 
children 

V Douet 1 , L Chang 1 , A Pritchett 1 , K Lee 1 , B Keating 1 , H Bartsch 2 , TL Jernigan 2 ' 3 , A Dale 2 ' 4 , N Akshoomoff 2 ' 3 , S Murray 5 , C Bloss 5 , 
DN Kennedy 6 , D Amaral 7 , J Gruen 8 , WE Kaufmann 9 , BJ Casey 10 , E Sowell 11 and T Ernst 1 for the Pediatric Imaging, Neurocognition and 
Genetics Study (PING) 

The neuregulin-1 {NRG1) gene is one of the best-validated risk genes for schizophrenia, and psychotic and bipolar disorders. The 
rs6994992 variant in the NRG1 promoter (SNP8NRG243177) is associated with altered frontal and temporal brain macrostructures 
and/or altered white matter density and integrity in schizophrenic adults, as well as healthy adults and neonates. However, the ages 
when these changes begin and whether neuroimaging phenotypes are associated with cognitive performance are not fully 
understood. Therefore, we investigated the association of the rs6994992 variant on developmental trajectories of brain macro- and 
microstructures, and their relationship with cognitive performance. A total of 972 healthy children aged 3-20 years had the 
genotype available for the NRG1-rs6994992 variant, and were evaluated with magnetic resonance imaging (MRI) and 
neuropsychological tests. Age-by-NRG1-rs6994992 interactions and genotype effects were assessed using a general additive model 
regression methodology, covaried for scanner type, socioeconomic status, sex and genetic ancestry factors. Compared with the 
C-carriers, children with the TT-risk-alleles had subtle microscopic and macroscopic changes in brain development that emerge or 
reverse during adolescence, a period when many psychiatric disorders are manifested. TT-children at late adolescence showed a 
lower age-dependent forniceal volume and lower fractional anisotropy; however, both measures were associated with better 
episodic memory performance. To our knowledge, we provide the first multimodal imaging evidence that genetic variation in NRG1 
is associated with age-related changes on brain development during typical childhood and adolescence, and delineated the altered 
patterns of development in multiple brain regions in children with the T-risk allele(s). 
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INTRODUCTION 

The neuregulin-1 gene {NRG1) encodes many structural and 
functional isoforms through alternative promoter usage (type 
l-IV). 1,2 These isoforms are involved in neuronal and glial 
development and synaptic plasticity. 3-6 Variants in NRG1 are 
associated with the diagnoses of schizophrenia and psychiatric 
disorders, 7-17 and contribute to impaired brain functions in those 
patients. 18 ' 19 The rs6994992 (SNP8NRG243177) is one of the most 
studied genes associated with increased risk for schizophrenic and 
psychotic symptoms, as well as cognitive deficits in 
schizophrenics. 20-22 This specific variant was also suggested to 
predict transition to psychosis in three independent cohorts with 
very high risk for psychosis. 20,23,24 Moreover, genetic variations at 
rs6994992 altered the promoter activity of the type IV NRG1 
isoform, which is expressed specifically in the brain and especially 
during the fetal period. 25 Higher transcript levels of type IV NRG1 
were also found in post-mortem brains of both healthy and 



schizophrenic adults with the TT-risk alleles. 2 Variation at 
rs6994992 in the type IV NRG1 results in differential expression 
that may have a critical role in the schizophrenia and psychosis 
pathology. 

Several magnetic resonance imaging (MRI) studies reported that 
NRG1-rs6994992 is associated with alterations in brain volumes 
and white matter integrity. In healthy adults, the T-risk allele is 
associated with smaller gray matter volumes in several gyri, 26 
higher white matter density in the right anterior internal capsule 27 
and lower white matter integrity in the left anterior thalamic 
radiation. 22,26 Furthermore, adults with first-onset schizophrenia 
had T-allele dose-dependent increases in lateral ventricle 
volumes. 28 Similar enlargements of the lateral ventricles along 
with altered developmental trajectories of gray and white matter 
were observed in TT-children (age 8-20 years) with childhood- 
onset schizophrenia compared with the C-carriers. 29 In contrast, 
healthy neonates with the rs6994992-TT and parental psychiatric 
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history (that is, bipolar disorder, depression, drug dependence/ 
abuse) showed greater gray matter volumes in the frontal and 
temporal lobes than the C-carriers. 30 

Genetic influences on brain development can change with 
age 31 and be modified by the environment. 32 Therefore, we 
propose that modulation of the NRG1 expression through genetic 
variation at the rs6694992 influences brain development and may 
show a pleiotropic effect in different brain structures, as well as an 
antagonistic pleiotropic effect at different age during normal brain 
development. To our knowledge, no study evaluated whether the 
functional rs6994992 variant influences age-related brain devel- 
opment from early childhood to late adolescence in a large 
population and at micro- and macrostructural levels. 

Cross-sectional studies across a large age span make it possible 
to evaluate age-dependent trajectories of brain development in 
relation to this particular genetic variation, and may improve our 
understanding of when these gene-related brain changes occur. 
The period of brain maturation in adolescence and early 
adulthood is a vulnerable time during which many mental 
illnesses are first expressed. Therefore, we investigated macro- 
and microstructural brain changes, using MRI, in a large cohort of 
healthy children aged 3-20 years. We focused on brain regions 
known to be associated with schizophrenia pathology, and 
reportedly affected by the rs6994992 variant (that is, mainly in 
the temporal, frontal and parietal lobes). Specifically, we 
hypothesized that compared with CC-alleles children, those with 
TT-alleles have differential age-dependent trajectories of micro- 
and macro-brain structures in these regions. 

Furthermore, a recent meta-analysis found that schizophrenic 
patients consistently showed slower processing speed and poorer 
episodic memory than healthy controls. 33 Therefore, we further 
evaluated whether the brain changes identified in the TT-children 
correlate with performance in these two domains. 

MATERIALS AND METHODS 

Data were obtained from the Pediatric Imaging Neurocognition and 
Genetic (PING) study database (http://ping.chd.ucsd.edu). PING is a cross- 
sectional normative data set from 1400 typically developing children aged 
3-20 years that was collected across ten academic institutions in the 
United States. The data set includes whole genome single-nucleotide 
polymorphism genotype, developmental and neuropsychological assess- 
ments (NIH Toolbox) and high-resolution brain MRI. The human research 
protection programs and institutional review boards of the participating 
institutions approved the design, experiments and consents (oral and 
written) used for the PING study. 

Subjects 

A total of 972 children from the PING database were included in this study. 
Details on recruitment and subject characteristics were reported 
previously. 34-37 Informed consent was obtained either from the partici- 
pants (^18 years old), or from their parents (<18 years of age) with child 
assent (7-17 years of age). Healthy children were screened to ensure that 
they had no self-reported or parent-reported history of major develop- 
mental, psychiatric or neurological disorders, brain injury or other medical 
conditions that would affect their brain development. Of the total, 250 
children were additionally examined by a physician to ensure their healthy 
status before enrollment. Individuals born prematurely (<36 weeks 
gestational age), with significant prenatal drug or alcohol exposure or 
with contraindications for MRI (for example, metallic or electronic implants, 
claustrophobia, dental braces or pregnancy) were also excluded. 

Genotyping, genotype imputation and genetic ancestry factor 
Genomic DNA was extracted from saliva for genotyping using the lllumina 
Human660W-Quad BeadChip. Replication and quality control filters (that is, 
sample call rate >99%, call rates >95%, minor allele frequency >5%) were 
performed 38 To assess reproducibility, about 1% of the samples were 
genotyped in duplicates to compare genotypes. The rs6994992 variant 
located on chromosome 8 in the promoter region of the neuregulin-1 
gene was imputed using MACH1 .0 39 with the default settings and a quality 
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control (R>0.5) for all individuals from the PING database as of 6 April 201 3. 
To control for race/ethnicity on the outcome measures, genetic ancestry 
factors (GAFs) were estimated using a standard method that includes 
multidimensional scaling analysis and actual estimates of local ancestry in 
admixed populations (LAMP) 40 as previously described 37 

Image acquisition and image processing 

The neuroimaging techniques, data acquisition and analyses were 
reported previously 34-36 and are available on the PING website (http:// 
ping.chd.ucsd.edu). Briefly, MRI was performed on ten 3 Tesla scanners 
from three manufacturers (Siemens, Erlangen, Germany; General Electric 
Medical, Milwaukee, Wl, USA; two Philips Medical, Andover, MA, USA). On 
the Siemens scanners, MRI included a three-dimensional T1 -weighted 
structural scan (TR/TE/TI = 21 70/2.78/1 100 ms, 7-degree flip angle, 
1. Ox 1.0 x 1.2 mm 3 resolution, 8 min acquisition time), a T2-weighted 
volume and a set of diffusion-weighted scans (dual spin-echo echo-planar 
imaging, TR/TE = 9500/91 s, 2 averages, 2.5 mm isotropic resolution, b- 
values = 1000s mm -2 , 30 diffusion directions, 10 min acquisition time). 
Imaging protocols for the other scanner platforms closely matched those 
for the Siemens scanners. Prospective motion compensation 41 was used 
for 74% of the T1 -weighted scans at UCSD 34 and 10.5% of those in Hawaii. 
After transferring the MRI data to the central processing site (UCSD) in 
DICOM format, all the images were automatically checked for protocol 
compliance on the basis of DICOM header information, and visually 
inspected for artifacts, including motion artifacts, signal inhomogeneities 
and ghosting. Scans rated as unacceptable on any rating were excluded 
from the study. 

Image analyses were performed using a modified Freesurfer software 
suite (http://surfer.nmr.mgh.harvard.edu/) that is based on routine three- 
dimensional T1 -weighted MRI volumes. FreeSurfer provides a method for 
automatic panellation of the cortical surface, by which each location on 
the surface is assigned a neuroanatomical label on the basis of 
probabilistic information estimated from a manually labeled training set. 
Details on surface-based cortical reconstruction and subcortical volumetric 
segmentation were reported previously 34 Delineation and labeling of the 
white matter tracts were performed with the atlas-based method, 
AtlasTrack, using diffusion tensor images (DTIs), T1 -weighted images and 
a probabilistic atlas of fiber tract locations and orientations 42 Validated for 
use in pediatric populations, AtlasTrack provides estimates for each tract 
the average mean diffusion, transverse and longitudinal diffusivity (TD and 
LD, respectively) and fractional anisotropy (FA). 

Fiber tract volumes were calculated as the number of voxels with fiber 
probability estimates greater than 0.08. 

Of the 972 children with rs6994992 imputed, all met quality criteria for 
cortical and subcortical measurements (thickness, area, volume), and 891 
for DTI measures (81 children had DTI scans that did not meet the quality 
criteria, and were excluded from our analysis. Out of the 81 children, 43 TC/ 
28 CC/10TT, Hardy-Weinberg equilibrium / = 0.48, P = 0.49). However, as 
we controlled for socioeconomic status, sex, device number and GAF in all 
the statistical models, only 897 children were included in the cortical and 
subcortical analysis and 822 for the DTI analysis. 

Episodic memory and processing speed 

The NIH Toolbox Cognition Battery evaluates six major cognitive 
domains 37,43 We analyzed data from two tasks that represent the two dom- 
ains implicated previously in schizophrenia: episodic memory (picture 
sequence memory test) and processing speed (pattern comparison process- 
ing speed test) 33 The picture sequence memory test first presents a 
sequence of thematically related pictures to the child on a computer screen, 
and then requires the child to sort a sequence of reshuffled pictures back to 
the order that was shown previously. Six to 18 items in each sequence over 
three trials were presented depending on the child's age (outcome: 
cumulative correct responses from a maximum of 18-48). For the pattern 
comparison processing speed task, the child has to determine whether two 
side-by-side pictures are similar by touching the smiley face or a sad face 
icon (ages 3-7 years), or a 'yes' or 'no' button (>7 years old), response on the 
screen (outcome: correct responses within the 90s time limit). 

Statistics 

Statistical analyses were performed using a customized user interface for 
the PING database, called 'Data Portal' (https://ping-dataportal.ucsd.edu/). 
The 'Data Portal' enables statistical analysis of genotype- and genotype-by- 
age effects on neural architectural phenotypes (from the imaging data), 
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Table 1. Participant characteristics (all values are in mean±s.e.) 



TT 



TC 



CC One-way ANOVA or 

P-value 



Children characteristics 

NRG1 genotype frequency, n (%) 
Age (years) 
Boys/girls 
GAF_Europe (%) 
GAF_Africa (%) 
GAF_American Indian (%) 
GAF_East Asia (%) 
GAF_Oceania (%) 
GAF_Central Asia (%) 



154 (15.84) 
11. 87 ±0.41 
80/74 
0.53 ±0.03 (29) 
0.90 ±0.02 (79) 
0.02 ±0.005 (17) 
0.30 ±0.03 (51) 
0.01 ±0.003 (38) 
0.05 ±0.01 (56) 



484 (49.79) 
12.11 ±0.23 
258/226 
0.63 ±0.02 (36) 
0.11 ±0.01 (10) 
0.04 ±0.005 (33) 
0.1 8 ±0.01 (30) 
0.008 ±0.001 (31) 
0.03 ±0.006 (33) 



334 (34.36) 
12.10 ±0.26 
171/163 
0.66 ±0.02 (35) 
0.1 3 ±0.01 (11) 
0.06 ±0.007 (50) 
0.1 2 ±0.01 (19) 
0.008 ±0.002 (31) 
0.01 ±0.005 (11) 



0.38 (HWE P-value) 
0.87 (F 0 . 14/ 971,2) 

0.83 (^.36, 972,2) b 

0.001 (F 6 . 68 , 970 , 2 ) a 

0.24 (F-i.42, 97o, 2 ) a 

0.0005 (F 7 . 61/ 970/2 ) a 

< 0.0001 (F 16 . 52 , 970,^ 
0.30 (FL2Q, 970,2)' 
0.06 (F 2 .91, 970 (2 ) 



Parent/guardian characteristics 

Household income (1 = < $5 K, 6 = 40 K-50 K, 12 = ^ $300 K) 6.78 ±0.1 9 6.84 ±0.11 6.62 ±0.1 4 0.44 (F 0 . 8 i, 928 , 2 ) a 

Highest education (7 = professional, 4= high school graduate, 5.91 ±0.09 5.91 ±0.05 5.64 ±0.06 0.002 (F 631 ' ; 938i2 Y 
1 = < 7 years of school) 

Highest occupation ( 7 = higher executives, 4 = clerical & sales 5.1 9 ± 0.1 3 5.1 8 ± 0.308 4.99 ± 0.09 0.23 [?^ A7i 89 2, 2 ) a 
workers, 1 = unskilled employee) 



MR scanners 

General Electric Medical 25 98 66 0.39 {& A[ 

Philips Medical 18 79 54 

Siemens 111 307 214 



Abbreviations: ANOVA, analysis of variance; HWE, Hardy-Weinberg equilibrium; MR, magnetic resonance. a F-value, sample size and degree of freedom. b x2 
value, sample size and degree of freedom. CC, TC and TT are genotypes. Significant genotype effects were observed only on the genetic ancestry factor (GAF) 
and highest education. 



and on relationships between neural phenotypes and neuropsychological 
measures using a general additive model 44 as implemented by the R 
program (http://www.r-project.org/). General additive model is a multiple 
linear regression model that allows smooth functions of the predictor 
variables or covariates where the degree of smoothness is data-driven 
rather than pre-specified 44 Each statistical model used age as a smooth 
independent variable, and included a linear term for the genotype and a 
smooth interaction between age and genotype. Also, all models accounted 
for group differences in socioeconomic status (reflecting the highest 
parental education and household income), sex, GAF 37 and study site. 
Cohen's f 2 effect sizes were calculated to estimate the contribution of 
smooth age-by-genotype interactions. Cortical area and thickness were 
additionally adjusted for the total mean 45 and subcortical volumes for 
intracranial volume. As cortical volumes are derived from regional areas 
and thickness, cortical volumes will not be reported to minimize multiple 
comparisons. For the vertex-based models, significance maps for cortical 
area and thickness were thresholded using the false discovery rate to 
correct for multiple comparisons at the 5% level. 46 For the region of 
interest (ROI)-based models, 44 cortical structures (22 for each hemisphere) 
were considered in this study, including mainly the temporal, frontal and 
parietal lobes. To correct for multiple comparisons we used Holm- 
Bonferroni sequential correction that involved a stepwise rejective multiple 
test procedure, testing the most extreme P-value against the strictest 
criterion, and the others against progressively less strict criteria 47,48 Results 
on the areas or thickness of the cortical structures were then considered 
significant when adjusted P-values were < 0.003 (a = 0.05). Five sub- 
cortical ROIs and lateral ventricles showed no hemispheric differences and 
were averaged before evaluating genotype and age-by-genotype effects. 
The subcortical ROI results were considered significant when adjusted 
P-values < 0.01 (a = 0.05). Likewise, as the measures for the six white 
matter tracts that were shown to be abnormal in patients with 
schizophrenia or altered in those with the T-risk allele (ALIC, ATR, fornix, 
corpus callosum, superior longitudinal fasciculus, uncinate) 20,27 showed 
no hemispheric differences, the values were averaged and adjusted 
P-values < 0.01 (a = 0.05) were considered significant. 

RESULTS 

Participant characteristics 

This study involved 972 children (mean age 12.07 ±4.97 years; 462 
girls and 510 boys; Table 1). The genotype distribution for the 



rs6994992 variant across all subjects did not deviate from the 
Hardy-Weinberg equilibrium (P = 0.38). Children in the three allele 
groups were similar in age, sex proportion, household income and 
parents/guardians' highest level of occupation. However, the 
genotype groups differed by GAFs and the highest level of 
education of the parents/guardians. As expected, the frequency of 
TT genotype was lower than TC or CC genotypes in children with 
European and American Indian GAF and higher in children with 
East-Asia or with Central-Asia GAF. Hence, GAF was controlled for 
in each statistical model. 

The rs6994992-T- allele is associated with altered development of 
subcortical structures 

Consistent with prior reports, 49,50 the volumes of pallidum, 
putamen, hippocampus and thalamus increased with age, but 
no genotype effects were found for these structures 
(Supplementary Table 1). On DTI, FA increased with age in basal 
ganglia structures (pallidum, caudate, putamen) and thalamus, but 
not in the hippocampus. Conversely, diffusivity (mean diffusion, 
LD and TD) generally showed age-related decreases in these 
subcortical structures. Age-by-genotype interactions were primar- 
ily found on FA in the basal ganglia. The age-related trajectories 
were significantly altered in the T-carriers for LD in the thalamus 
(Figure 1) and for TD in the putamen (Supplementary Figure 1) 
and pallidum (Figure 1). Furthermore, T-carriers had a higher FA in 
the pallidum than CC-carriers during adolescence but lower FA 
during younger childhood (Figure 1). Similar patterns were also 
found for the caudate and putamen (Supplementary Figure 1). 
Conversely, TD in these regions was higher in T-carriers than CC- 
carriers at younger ages, but showed a lower age-related TD 
during adolescence (Figure 1). Similar to schizophrenic adults, 28 
healthy TT-adolescents (>14 years of age) showed a trend for 
larger and steeper age-dependent increase in the lateral ventricles 
compared with C-carriers, leading to 37% larger estimated 
volumes at 20 years of age (Figure 1). Our statistical model 
explained 23.4% of the variance of the lateral ventricle volumes 
(age-by-genotype interaction, P = 0.02, R 2 (adjusted) = 0.21, effect 
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Figure 1. Age-related differences in DTI metrics of subcortical structures (top row) and in lateral ventricle volume (bottom row). Percentages 
are expressed as increases (sign +) or decreases (sign -) in values between the TT-children (green) and the CC-children (blue). Black vertical 
arrows indicate when the age-by-genotype interactions emerge or reverse. The bottom right panel illustrates the location of subcortical 
structures of interest, as well as a statistical P-value map for FA generated with the ROI-based model. DTI, diffusion tensor image; FA, fractional 
anisotropy; LD, longitudinal diffusivity; ROI, region of interest; TD, transverse diffusivity. 



size (/ 2 ) =0.266, power >0.90, a = 0.05; a model excluding 
genotype explained only 11% of the variance). 

The rs6994992-T-risk allele alters development of cortical 
morphometry 

Figure 2 summarizes age-by-genotype interactions and genotype 
effects on developmental trajectories of the cortex using both 
vertex-based (Figure 2a) and ROI-based (Figure 2b) models. On 
both models, age-by-genotype interactions were found predomi- 
nantly in the temporal, parietal and frontal lobes. For instance, the 
vertex-based model showed age-by-genotype interactions in 
thickness for the right supramarginal gyrus, the right superior 
temporal gyrus and the left fusiform, and in area for the middle 
frontal, inferior parietal right caudal anterior cingulate and left 
precuneus (Figure 2a). Most of the cortical phenotypes affected by 
the NRG1 variant emerged or reversed during adolescence 
(Supplementary Figures 2-5). 

Three general patterns of age-by-genotype interaction were 
observed with the ROI-based models (Figure 3). The first pattern 
is exemplified by the rostral middle frontal gyrus surface area. 
Young TT-children did not exhibit genotype differences, but TT- 
adolescents showed a lesser age-related decline in area (Figure 3a) 
relative to C-carriers. This pattern was also observed for the area of 
the right frontal pole (Supplementary Figure 2). 

The second pattern involved a reversal of neuroimaging 
phenotypes with age, commonly with a crossing of trajectories 
during adolescence. For instance, in the right caudal middle 
frontal gyrus, young TT-children had a larger area than C-carriers, 
whereas TT-adolescents had a smaller area (Figure 3b). Similar 
interactions were found for the area of the lateral orbitofrontal 



cortex and the thickness of the left supramarginal gyrus 
(Supplementary Figures 3). Conversely, in the left caudal anterior 
cingulate, younger TT-children had smaller area (Figure 3c) and 
thickness (Supplementary Figure 3) than C-carriers, whereas the 
differences were reversed in adolescents. The right superior 
frontal gyrus exhibited this pattern on thickness, and the right 
inferior temporal gyrus on area (Supplementary Figure 3). 

The third pattern was found for the right superior temporal 
gyrus, where the age-related curve of the TT-carriers for area 
followed a U-shape with a trough at 12.5 years of age, whereas a 
slightly inverted U-shape was observed in the C-carriers 
(Figure 3d). This pattern was also found for the thickness of the 
right inferior parietal cortex (Supplementary Figure 4). 

Independent of age and compared with CC-carriers, TT-children 
had larger area of the right inferior parietal cortex (Supplementary 
Figure 5) and showed a steeper decrease in area of the left 
precuneus (Supplementary Figure 5). 

The rs6994992-T-risk allele is associated with alterations in 
forniceal white matter 

Among the white matter tracts evaluated (ALIC, ATR, Fornix, 
Corpus callosum, SLF), an age-by-genotype interaction was found 
only in the fornix volume and DTI metrics (Figure 4). TT-children 
showed inverted U-shapes in the growth of the volume (peaking 
at 11.5 years) and FA (peaking at 14 years) of the fornix. 
Conversely, CC-carriers showed continued growth or a plateau in 
both metrics during their second half of adolescence (Figures 4a 
and b). These differential developmental trajectories resulted in 
smaller volume (-6.1%) and lower FA (-3.4%) in TT-children 
compared with CC-children by age 20. 



Translational Psychiatry (2014), 1 -9 



© 2014 Macmillan Publishers Limited 



Effect of rs6994992 variant on brain development 
V Douet et al 



P-maps for age x NRG1 interactions 
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Figure 2. Summary of NRG1 genotype main effects and age-by-genotype interactions in cortical measures, (a) Statistical P-value maps 
showing significant interactions between smooth age and the three NRG1 genotypes (TT f TC and CC) in healthy children. The maps were 
generated from general additive models with cortical area and thickness as the dependent variable at each location (vertex) across the 
surface. Sex, device, GAF, SES were included as covariates, and findings were corrected for multiple comparisons with false discovery rate 
(FDR; a = 0.05). The FDR threshold was obtained for left and right hemispheres combined and based on the cortical structures maps, (b) Table 
of NRG1-rs6994992 genotype differences and genotype-by-age interactions on cortical structures. The P-values (two-way analysis 
of covariance) of regions showing statistical differences with both the ROI-based and the vertex-based models are indicated in italic. The 
P-values that survived Holm-Bonferroni correction are in bold. GAF, genetic ancestry factor; ROI, region of interest; SES, socioeconomic status. 



Paradoxically, these TT-children with smaller volumes and lower 
FA of the fornix actually had higher picture sequence memory test 
scores (better episodic memory; Figure 4c). No significant 
correlations were found between pattern comparison processing 
speed test scores (processing speed) and brain imaging measures 
in the fornix. 



DISCUSSION 

To our knowledge, this is the first demonstration that the 
functional genetic variant rs6994992 in the NRG1 gene is 
associated with both microscopic and macroscopic alterations in 
brain development of typically developing healthy children and 
adolescents. Although the changes are subtle, some of our 
findings on cortical structures are consistent with those reported 
in neonates with high familial risk for psychosis 30 and in 
schizophrenic adults and adolescents. 26-28 This is also the first 
report to show that most of the NRG1-rs6994992 effects on 



neuroimaging phenotypes emerge or reverse during adolescence. 
Therefore, the neuroimaging phenotypes associated with this risk 
gene appear to evolve coincidentally with mental disorders that 
often emerge during adolescence, most likely due to the 
differential expression of the type IV NRG1 mediated by variation 
at rs6994992. 

The NRG1-rs6994992 variant is located in the alternative 
promoter region of type IV NRG1, and those with the TT-alleles 
produce higher transcript levels (up to 40% higher) of the type IV 
NRG1 than those with the CC-alleles, as shown in the post-mortem 
brains of both healthy and schizophrenic individuals 2 Likewise, 
the NRG1 promoter of the T-alleles is 65% more active than that of 
C-alleles in cell cultures 25 Therefore, rs6994992-TT children might 
produce more type IV NRG1 during brain development, which 
could lead to structural changes that ultimately increase the risk 
for mental disorders. However, the associated cellular mediators 
on structural MRI measures are unknown, and thus the molecular 
mechanism implicating type IV NRG1 remains speculative. 
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Figure 3. Examples illustrating the different patterns of the age-by-genotype interactions on cortical area (obtained from the ROI-based 
model). The graphs exemplify age-related effects of NRG1 that emerge (a) or reverse (b and c) at adolescence, and a more complex age-by 
genotype interaction (d). Percentages are expressed as increases (sign +) or decreases (sign -) in values between the TT-children (green) and 
the CC-children (blue). Black vertical arrows indicate when the age-by-genotype interactions emerge or reverse. ROI, region of interest. 



Nevertheless, assessing how this variant affects brain develop- 
ment and whether the expression of type IV NRG1 through this 
variation is involved in this process may lead to new therapeutic 
targets. 

On cortical structures, both vertex- and ROI-based analyses 
showed that genetic variation at rs6994992 is associated primarily 
with alterations in the growth trajectories of the cortical area, and 
hence volume, in the frontal, temporal and parietal lobes. Our 
findings are consistent with a prior study that showed TT-allele 
children with childhood-onset schizophrenia had greater total 
gray matter volume and changes in shapes and forms of 
the frontal, temporal and parietal lobes. 29 Similar to our study, 
TT-neonates with high familial risk for psychosis had larger 
volumes than C-carriers primarily in the frontal and temporal 
lobes, and smaller volumes in several smaller regions across the 
brain. 30 In addition, our study demonstrates that variation at 
rs6994992 affects regions that are often associated with volume 
loss in adolescents with ultra-high risk for psychosis, such as that 
in the left caudal anterior cingulate 51 and left precuneus. 52 
Replicating findings of rs6994992 effects on cortical structures in 
such a large cohort of typically developing children strongly 
supports the biological relevance of the gene's effects on newly 
identified regions involved in psychosis, such as the right inferior 
parietal gyrus or the supramarginal gyrus. However, these prior 
studies did not investigate age-by-NRG1 genotype interactions, 
perhaps due to the smaller sample size and narrower age ranges. 

Larger ventricular volume in TT-carriers compared with 
C-carriers started in adolescence. Although the statistical sig- 
nificance for this trait is borderline (P = 0.02), similar ventricular 



enlargement associated with the rs6994992 variant was reported 
in schizophrenic adults and unaffected relatives, 28 as well as in 
children with childhood-onset schizophrenia (8-20 years of age) 29 
On subcortical structures, we did not find smaller basal ganglia 
and thalamus volumes as reported in TT-neonates. 30 However, 
A//?G7-related alterations were observed on striatal FA, an indicator 
of white matter coherence and axonal organization, 53 and on TD, 
which may reflect myelination or glial cell morphology. 54-56 
Therefore, the rs6994992 variant may affect white matter through 
myelination rather than axonal integrity, since the LD was 
normal. 55,57 In contrast, the rs6994992 variant may affect the 
axonal integrity or extracellular space of the thalamus, as indicated 
by changes in both the LD and FA. Again, these subcortical 
changes typically emerged at adolescence. These results highlight 
the pleiotropic role of the rs6994992-A//?G 7 T-risk allele on brain 
development, and suggest that peripubertal factors, such as 
hormones or environment, might further impact the gene 
expression. 

The developmental trajectories of FA, mean diffusion and 
volume of the fornix in the PING cohort resemble those in another 
smaller healthy cohort (5-83 years of age) that showed inverted 
U-shaped curves on forniceal FA and volumes, peaking near age 
20 years. 58 The magnitude of declines in forniceal volume and FA 
in our TT-children is also comparable with the findings in that 
cohort, 58 suggesting an accelerated aging effect of NRG1- 
rs6994992 T-allele on fornix metrics. 

As reduced fornix volume was found in both schizophrenic 
adolescents 59 and adults, 60 our finding of smaller fornix in TT-risk 
children suggests a presymptomatic endophenotype. The 
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Figure 4. Associations between rs6994992-NRG1 genotype and fornix measures over age. The fornix tract showed age and NRG1 genotype 
dependent trajectories in volume (a) and FA (b). (c) Correlation between episodic memory task performance and forniceal metrics as a 
function of the NRG1 genotype. The center panel illustrates the location of the fornix. Percentages in the graphs are expressed as increases 
(sign +) or decreases (sign -) in values between the TT-children (green) and the CC-children (blue). Black vertical arrows indicate when the 
age-by-genotype interactions emerge or reverse. The P-values that survived Holm-Bonferroni correction are in bold. FA, fractional anisotropy; 
PSMT picture sequence memory test. 



emergence of these alterations during adolescence suggests a 
possible association with hormonal changes that occur during 
puberty. However, longitudinal follow-up studies are required to 
assess whether these children will develop psychiatric symptoms 
or what factors might trigger the onset of symptoms. 

Finally, smaller fornix volume in our TT-children was associated 
with better episodic memory (picture sequence memory test), 
which is opposite to the poorer episodic memory found in 
schizophrenics. 33 In contrast, our findings suggest a positive 
impact of the rs6994992 variant on neuropsychological perfor- 
mance and are consistent with those reported by Keri. 61 The 
author described that the TT-carriers had the highest creative- 
thinking scores in healthy adults with high intellectual and 
academic performance. Creative thinking and episodic memory 
were previously reported to be related to psychosis. 20,23,24 The 
better performance in the healthy subjects may reflect an 
antagonistic pleiotropy effect, with better performance at younger 
age and poorer performance with aging, in these individuals with 
the T-risk allele. Nevertheless, abnormalities in fornix volume and 
diffusivity often precede the clinical deterioration of cognitive 
function and episodic memory in schizophrenics. 62-64 Therefore, 



our results further support the hypothesis that NRG1-TT alleles 
confer risk for the schizophrenia phenotype via lesser myelination 
and atrophy of the white matter, especially the fornix. 

Limitations 

A potential limitation with a large scale study is the wide diversity 
in genetic backgrounds of the participants. However, a major 
strength of the PING data set is that the study was designed to 
detect a genetic variant effect size of 0.1 with a power >0.99 and a 
sample size of at least 800 children, assuming a total of seven 
independent variables and a type I error rate of 1x10 -7 . 
However, findings from the PING study are somewhat limited by 
the cross-sectional study design. For instance, outcome measures 
reflect changes in neurophenotypes of the children at the time of 
the scan, although the phenotype could have been present earlier. 
Future longitudinal studies are needed to validate the role of the 
type IV NRG1 in neurodevelopmental trajectories, aging and 
pathology. To minimize erroneous conclusions about the effects 
of age and the NRG1-rs6994992rs6994992 variant on changes in 
brain structures over time, all analyses included fixed traits such as 
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sex, ethnicity and socioeconomic status as covariates. Another 
potential limitation is the possible influence of subclinical disease 
on brain measures in some of these children, despite our efforts to 
include only healthy children by careful screening and evaluations. 
As schizophrenia and psychosis may develop subclinical^ years 
before clinical symptoms are detectable, long-term follow-up 
examinations are needed to determine whether some of these 
individuals will develop these clinical disorders. 

CONCLUSION 

The rs6994992-polymorphism in the NRG1 gene is associated with 
divergent age-related changes in gray and white matter regions. 
Those changes are expressed during adolescence, a period when 
symptoms for neuropsychiatric disorder commonly arise. Our 
study bridges the existing gaps and validates prior smaller studies 
on the effects of rs6994992-/W?G 7 on developmental trajectories 
during childhood and adolescence, and ties in with prior findings 
from neonates through adulthood. Future research on the 
molecular mechanism underlying the roles of NRG1 and the 
genes involved in its pathway, such as ERBB4, on brain develop- 
ment and cognitive performance may ultimately lead to early 
detection and improved treatment of schizophrenia and psychotic 
disorders. 
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